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Present cosmological observations yield an upper bound on the neutrino mass 
which is significantly stronger than laboratory bounds. However, the exact value 
of the cosmological bound is model dependent and therefore less robust. Here, I 
review the current status of cosmological neutrino mass bounds and also discuss 
implications for sterile neutrinos and LSND in particular. 



1. Introduction 

The absolute value of neutrino masses are very difficult to measure experimen- 
tally. On the other hand, mass differences between neutrino mass eigenstates, 
(toi, m2, fna), can be measured in neutrino oscillation experiments. 

The combination of all currently available data suggests two important mass 
differences in the neutrino mass hierarchy. The solar mass difference o f Srn^ r, ~ 
7 x eV 2 and the atmospheric mass difference Sm^ — 2.6 x 10~ 3 eV 2 l ' (see 
also the contribution by C. Giunti to the present volume). 

In the simplest case where neutrino masses are hierarchical these results sug- 
gest that mi ~ 0, rn,o ~ Sm so i ar , and ~ ^^atmospheric- If the hierar- 
chy is inverted l 4 l v l v l ( l g l t' l one instead finds 7713 ~ 0, m-x ~ <5m a t mosp heric, ancl 
mi ~ 5m p Am nanh srj <•. ■ However , it is also possible that neutrino masses are degener- 
ate UJjyjJJ^UjiL^LjiLii^^ mi ~ rri2 ~ m?, <Sm a t mosp hcric, in which case 
oscillation experiments are not useful for determining the absolute mass scale. 

Experiments which rely on kinematical effects of the neutrino mass offer the 
strongest probe of this overall mass scale. Tritium decay measurements have been 
able to put an upper limit on the electron neutrino mass of 2.2-2.3 eV (95% conf.) I2H 
(see also the contribution by C. Kraus in the present volume). However, cosmology 
at present yields an even stronger limit which is also based on the kinematics of 
neutrino mass. 

Neutrinos decouple at a temperature of 1-2 MeV in the early universe, shortly 
before electron-positron annihilation. Therefore their temperature is lower than the 
photon temperature by a factor (4/11) 1 / 3 . This again means that the total neutrino 
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number density is related to the photon number density by 

9 

Massive neutrinos with masses m ^> T ~ 2.4 x 10~ 4 eV are no n-relativis tic at 
present and therefore contribute to the cosmological matter density E2ESE1 

" ~ 92.5 eV 1 ' 

calculated for a present day photon temperature T = 2.728K. Here, Yl m v = 
mi + mi + 777-3 ■ However, because they are so light these neutrinos free stream on 
a scale of roughly k ~ 0.03777 c v^m 2 h Mpc -1 ESEHEI Below this scale neutrino 
perturbations are completely erased and t here fore the matter power spectrum is 
suppressed, roughly by AP/P ~ -8fi„/fi m E3. 

This power spectrum suppression allows for a determination of the neutrino mass 
from measurements of the matter power spectrum on large scales. This matter 
spectrum is related to the galaxy correlation spectrum measured in large scale 
structure (LSS) surveys via the bias para meter , b 2 = P g (k)/P m (k). Such analyses 
have been performed several times before EHEHI most recently using data from the 
2dF galaxy survey ^4 

However, using large scale structure data alone does not allow for a precise 
determination of neutrino masses, because the power spectrum suppression can 
also be caused by changes in other parameters, such as the matter density or the 
Hubble parameter. 

Therefore it is necessary to add information on other para meters from the cosmic 
microwave background (CMB). This has been done in using ealier 

CMB d ata. More re cently the precise data from WMAP «™ has been used for this 
purpose E2E3E32SI to derive a limit of 0.7-1.0 eV for the sum of neutrino masses. 



2. Cosmological data and likelihood analysis 

The extraction of cosmological parameters from cosmological data is a difficult pro- 
cess since for both CMB and LSS the power spectra depend on a plethora of differ- 
ent parameters. Furthermore, since the CMB and matter power spectra depend on 
many different parameters one might worry that an analysis which is too restricted 
in parameter space could give spuriously strong limits on a given parameter. 

The most recent cosmological data is in excellent agreement with a flat ACDM 
model, the only non-standard feature being the apparently very high optical depth 
to reionization. Therefore the natural benchmark against which non-standard neu- 
trino physics can be tested is a model with the following free parameters: f2 TO , the 
matter density, the curvature parameter, fib, the baryon density, Hq, the Hubble 
parameter, n s , the scalar spectral index of the primordial fluctuation spectrum, r, 
the optical depth to reionization, Q, the normalization of the CMB power spec- 
trum, 6, the bias parameter, and finally the two parameters related to neutrino 
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physics, VL^h? and N v . The analysis can be restricted to geometrically flat models, 
i.e. Q = Q m + CI \ = L For the purpose of actual power spectrum calculations, the 
CMBFAST package can be used. 



2.1. LSS data 

At present there ar e two l arge galaxy surveys of comparable size, the Sloan Digital 
Sky Survey (SDSS)l^l^and the 2dFGRS (2 degree Field Galaxy Redshift Survey) 
™4 Once the SDSS is completed in 2005 it will be significantly larger and more 
accurate than the 2dFGRS. At present the two surveys are, however, comparable 
in precision and here we discuss constraints from the 2dFGRS alone. 

Tegmark, Hamilton and Xul^have calculated a power spectrum, P(k), from 
this data, which we use in the present work. The 2dFGRS data extends to very 
small scales where there are large effects of non-linearity. Since we calculate only 
linear power spectra, we follow standard procedures and use only data on scales 
larger than k = 0.2 h Mpc -1 , where effects of non-linearity should be minimal (see 
for instance Ref. ^El for a discussion) . With this cut the number of data points for 
the power spectrum reduces to f 8. 



2.2. CMB data 

The CMB temperature fluctuations are conveniently described in terms of the spher- 
ical harmonics power spectrum C/ = (|et/ m | 2 ), where ^-(6,4>) — ^2 lm ai m Yi m (6 , (j>) . 
Since Thomson scattering polarizes light there are also power spectra coming from 
the polarization. The polarization can be divided into a curl-free (E) and a curl 
(B) component, yielding four independent power spectra: Ct,i,Ce,i,Cb,i and the 
temperature iS-polarization cross-correlation Cte,i- 

Th e WMAP exper iment have reported data only on Ct.i and Cte I, as described 



r exper 
4 45 46 



in Ref. EMM 

We have performed the likelihood analysis using the prescription given by 
ihe WMAP c ollaboration which includes the correlation between different Ci's 
| 33 | 43 | 44 | 45 | 46 | Foreground contamination has already been subtracted from their 
published data. 

In parts of the data analysis we also add other CMB data from the compilation 
by Wang et al. which includes data at high I. Altogether this data set has 28 
data points. 



3. Neutrino mass bounds 

The analysis presented here was originally published in Ref. an( j more details 
can be found there. 

We have calculated \ 2 as a function of neutrino mass while marginalizing over 
all other cosmological parameters. This has been done using the data sets described 
above. In the first case we have calculated the constraint using the WMAP Ct,i 
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combined with the 2dFGRS data, and in the second case we have added the polar- 
ization measurement from WMAP. Finally we have added the additional constraint 
from the HST key project and the Supernova Cosmology Project. It should also be 
noted that when constraining the neutrino mass it has in all cases been assumed 
that N v is equal to the standard model value of 3.04. Later we relax this condition 
in order to study the LSND bound. 

The result is shown in Fig. 1. As can be seen from the figure the 95% confidence 
upper limit on the sum of neutrino masses is ^2 l m u < 1.01 eV (95% conf.) using 
the case with priors. This value is completely consistent with what is found in 
Ref. ESI where simple Gaussian priors from WMAP were added to the 2dFGRS 
data analysis. For the three cases studied the upper limits on m v can be found 
in Table 1. 



Case £m„ (95% C.L.) 

1 1.01 eV 

2 1.20 eV 

3 2.12 eV 



In the middle panel of Fig. 1 we show the best fit value of Hq for a given Vl u h 2 . It 
is clear that an increasing value of Y2,m v can be compensated by a decrease in H . 
Even though the data yields a strong constraint on fl m h 2 there is no independent 
constraint on Q m in itself. Therefore, an decreasing Hq leads to an increasing Q m . 
This can be seen in the bottom panel of Fig. 1. 

When the HST prior on Hq is relaxed a higher value of m v is allowed, in the 
case with only WMAP and 2dFGRS data the upper bound is Vt v h 2 < 0.023 (95% 
conf.), corresponding to a neutrino mass of 0.71 eV for each of the three neutrinos. 

16 HI 

This effect was also found by Elgar0y and Lahav in their analysis of the 
effects of priors on the determination of ^2 m v- 

However, as can also be seen from the figure, the addition of high-Z CMB data 
from the Want et al. compilation also shrinks the allowed range of m v signif- 
icantly. The reason is that there is a significant overlap of the scales probed by 
high-Z CMB experiments and the 2dFGRS survey. Therefore, even though we use 
bias as a free fitting parameter, it is strongly constrained by the fact that the CMB 
and 2dFGRS data essentially cover much of the s ame range in fc-space. 

It should be noted that Elgar0y and Lahav ^3 fi n d that bias does not play 
any role in determining the bound on Tn v . At first this seems to contradict the 
discussion here, and also wha t was found from a Fisher matrix analysis in Ref. 
The reason is that in Ref. redshift distortions are included in the 2dFGRS data 
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analysis. Given a constraint on the amplitude of fluctuations from CMB data, and 
a constraint on £l m h 2 , this effectively constrains the bias parameter. Therefore 
adding a further constraint on bias in their analysis does not change the results. 

Neutrinoless double beta decay - Recently it was claimed that the Heidelberg- 
Moscow experiment yields positive evidence for neutrinoless double beta decay. 
Such experiments probe the 'effective electron neutrino mass m ee = | U^m v .\. 
Given the uncertainties in the involved nuclear matrix elements the Heidelberg- 
Moscow result leads to a mass of m ee = 0.3 — 1.4 eV. If this is true then the mass 
eigenstates are necessarily degenerate, and ~ 3m ee . Taking the WMAP 

result of m v < 0-70 e V at face value seems to be inconsistent with the Heidelberg- 
Moscow resulfESl. However, already if Ly-a forest data and a constraint on the bias 
parameter is not used in the analysis the upper bound of J2 m f — 1-01 e V is still 
consistent. For this reason it is probably premature to rule out the claimed evidence 
for neutrinoless double beta decay. 

Evidence for a non-zero neutrino mass - In a recent paper it was noted that 
there is a preference for a non-zero neutrino mass if a measurement of the bias 
parameter from X-ray clusters is added to the CMB and large scale structure data. 
This result arises because the X-ray data prefers a low value of as (bias), which 
is incompatible with the WMAP and 2dF result at the 2a level. While this is an 
interesting finding it is clear that the X-ray data is subject to a serious problem with 
systematic uncertainties, such as the calibration of the mass-temperature relation. 
Therefore the result more likely points to a problem with the interpretation of the 
X-ray data than to evidence of a non-zero neutrino mass. 

4. Sterile neutrinos 

In Ref. El it was shown that there is a degeneracy between the neutrino mass 
{^2,rn v ) and the relativistic energy density, parameterized in terms of the effective 
number of neutrino species, N v . 

As can be seen from Fig. 2, the best fit actually is actually shifted to higher 
^2 m v when N v increases, and the conclusion is that a model with high neutrino 
mass and additional relativistic energy density can provide acceptable fits to the 
data. As a function of N v the upper bound on m v (at 95% confidence) can be 
seen in Table 2. 



effective N v 


£m„ (95% C.L.) 


3 


1.01 eV 


4 


1.38 eV 


5 


2.12 eV 
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Figure 1. The top panel shows x 2 as a function of Yl m v for different choices of priors. The 
dotted line is for WMAP + 2dFGRS data alone, the dashed line is with the additional Wang et al. 
data. The full line is for additional HST and SNI-a priors as discussed in the text. The horizontal 
lines show A\ 2 = 1 and 4 respectively. The middle panel shows the best fit values of Ho for a 
given J] m v . The horizontal lines show the HST key project lcr limit of Ho = 72±8 kmsMpc . 
Finally, the lower panel shows best fit values of Q m - In this case the horizontal line corresponds 
to the SNI-a lcr upper limit of O m < 0.42. 



This has significant implications for attempts to constrain the LSND experi- 
ment using the present cosmological data. Pierce and Murayama concl ude from the 
present MAP limit that the LSND result is excluded^ (see also Ref.^Hj). 

However, for several reasons this conclusion does not follow trivially from 
the present data. In general the three mass differences implied by Solar, at- 
mospheric and the LSND neutrino m eas urements can be arranged into either 
2+2 or 3+1 schemes. Recent analyses 1^ of experimental data have shown that 
the 2+2 models are ruled out. The 3+1 scheme with a single massive state, 
ni4, which makes up the LSND mass gap, is still marginally allowed in a few 
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small windows in the (Am 2 , sin 2 29) plane. These gaps are at (Am 2 , sin 2 20) ~ 
(0.8eV 2 ,2 x lCT 3 ),(1.8eV 2 ,8 x 10~ 4 ), (6eV 2 , 1.5 x lO" 3 ) and (10eV 2 ,1.5 x lO" 3 ). 
These four windows corresponds to masses of 0.9,1.4,2.5 and 3.2 eV respectively. 
From the Solar and atmospheric neutrino results the three light mass eigenstates 
contribute only about 0.1 eV of mass if they are hierarchical. This means that the 
sum of all mass eigenstate is close to m^. 

The limit for N v — 4 which corresponds roughly to the LSND scenario is ^2 m v < 
1.4 eV, which still leaves the lowest of the remaining windows. The second window 
at m ~ 1.8 eV is disfavoured by the data, but not at very high significance. 




Figure 2. Ax 2 as a function of m w for various different values of N v . The full line is for 
N„ = 3, the dotted for N„ = 4, and the dashed for N v = 5. A% 2 is calculated relative to the best 
fit N v = 3 model. 



5. Discussion 

We have calculated improved constraints on neutrino masses and the cosmological 
relativistic energy density, using the new WMAP data together with data from the 
2dFGRS galaxy survey. 

Using CMB and LSS data together with a prior from the HST key project on Hq 
yielded an upper bound of 'Y^m v < 1.01 eV (95% conf.). While this excludes most 
of the parameter range suggested by the claimed evidence for neutrinoless double 
beta decay in the Heidelberg-Moscow experiment, it seems premature to rule out 
this claim based on cosmological observations. 
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Another issue where the cosmological upper bound on neutrino masses is very 
important is for the prospects of directly measuring neutrino masses in tritium end- 
point measurements. The successor to the Mainz experiment, KATRIN, is designed 
to measure an electron neutrino mass of roughly 0.2 eV, or in terms of the sum of 
neutrino mass eigenstates, ^2m u < 0.75 eV (see contribution by Guido Drexlin to 
the present volume). The WMAP result of J2 m v — 0-7 eV (95% conf.) already 
seems to exclude a positive measurement of mass in KATRIN. However, this very 
tight limit depends on priors, as well as Ly-a forest data, and the more conservative 
present limit of < 1.01 eV (95% conf.) does not exclude that KATRIN will 

detect a neutrino mass. 

Finally, we also found that the neutrino mass bound depends on the total number 
of light neutrino species. In scenarios with sterile neutrinos this is an important 
factor. For instance in 3+1 models the mass bound increases from 1.0 eV to 1.4 eV, 
meaning that the LSND result is not ruled out by cosmological observations yet. 
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